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Abstract 

We investigate relations between spectral properties of a single-centre point-interaction 
Hamiltonian describing a particle confined to a bounded domain Q, C R d , d — 2,3, with 
Dirichlet boundary, and the geometry of Q. For this class of operators Krein's formula yields 
an explicit representation of the resolvent in terms of the integral kernel of the unperturbed 
one, (— + z) 1 . We use a moving plane analysis to characterize the behaviour of the 
ground-state energy of the Hamiltonian with respect to the point-interaction position and the 
shape of Q, in particular, we establish some conditions showing how to place the interaction 
to optimize the principal eigenvalue. 

1 Introduction 

Relations between geometry of a domain and spectral properties of corresponding operators 
belong to the most traditional question in mathematical physics; one can recall, e.g., the Faber- 
Krahn inequality [Fa23, Kr25 or the Payne-Polya- Weinberger conjecture [PPW55] proved by 
Ashbaugh and Benguria [AB92a, AB92b . A more recent example concerns the situation where 
the domain in question is not simply connected and one asks, in particular, how to place a circular 
hard-wall obstacle within a circular planar cavity to minimize the ground-state eigenvalue; it 
appears that the minimum is reached when the obstacle touches the boundary [HKKOlj . 

In connection with the last mentioned problem one can also ask what happens if such a 
"hard" obstacle is replaced by another object, say, by a potential barrier or well. In this paper 
we are going to address this question in the particular case when such a potential is singular, 
in other words, a point interaction. Recalling basic results about these interactions AGHH05 
we see that the problem makes sense in dimension d < 3. Furthermore, a simple perturbative 
argument shows that in the one dimensional situation the answer may depend on the sign of the 
(5 potential; we restrict here our attention to the more singular case of dimension d = 2, 3. 

Problems of this type were to our knowledge solved so far only in cases where the domain 



has a simple geometry such as a straight strip in K , see [EGST96], or a planar layer in R , see 
EN02], where the eigenvalue problem can be solved more or less explicitly. Here we consider 
that the domain f2, to which the particle is confined, is bounded and otherwise quite general, 
see below. The operator of interest will be the corresponding Dirichlet Laplacian perturbed by 
a single point interaction of a fixed coupling constant; we will ask about the dependence of its 
principal eigenvalue on the perturbation position. Using a figurative expression, to be made 
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precise below, we are going to show that the ground-state energy increases as the point moves 
towards the boundary of SI, in contrast to the case of a Dirichlet obstacle mentioned above. 

Our method is based on the fact that one is able to express the resolvent of the operator in 
question by means of the Krein's formula. The resulting spectral condition allow us to character- 
ize the principal eigenvalue as a function of the interaction position in a non perturbative setting. 
Then we exploit the maximum principle and a domain reflection technique, analogous to the one 
used in [HKKOlJ, to demonstrate our main result, Theorem 14.31 expressing strict monotonicity 
of the principal eigenvalue with respect to certain directions. Using this conclusion, we formulate 
then some conditions under which the eigenvalue reaches its minimum value for a fixed SI. 



2 Confined point interactions in two and three dimensions 

The definition domain and the spectral properties of point interaction Hamiltonians in dimensions 
two and three are usually expressed in terms of the free Green function, i.e. the integral kernel 
of the operator (—A + z) , cf. [AGHHQ5] . This is true both if the configuration space is the 
whole M. d or if the particle is confined to a subset SI of it by a hard wall corresponding to Dirichlet 
boundary condition. In the latter case Green's function is defined by the equation 

r (-A + z)gz(x,x!) = s(x- x!) 

< ; x' e Q. (2.1) 

which admits a solution in L 2 (fl) whenever — z does not belong to the spectrum of the Dirichlet 
Laplacian — A^ 1 defined in the standard way through the associated quadratic form [RS78J. 
Throughout this paper, SI is supposed to be an open regular set in M d , d = 2, 3, bounded and 
connected, not necessarily simply; we assume that dtt is piecewise C 1 . 

Under these hypotheses — A^ has a purely discrete spectrum; we denote by {A„} ngNo , where 
No := NU{0}, its eigenvalues, and by {VVfcl^i the corresponding system of eigenfunctions, 

nGNo 

N n being the multiplicity of the n-th eigenvalue. Projecting (|2.1|) on the vectors ip n ,k, it is easy 
to check the validity of the following standard Fourier expansion of G§ , 

g^)= V t^ifMEl. (2 . 2) 

nGN ™ 
k<N n 

An alternative representation of the kernel, which will be extensively used in this paper, can 
be given in terms of the 'free' Green's function G z (x,af), in other words, the integral kernel of 
(— A + zy 1 in the whole spac^E 

g z (x,x) = ^-K (^\x-x!\) inM 2 (2.3) 
Zir 

-y/z\x-x'\ 

Q z {^ll) = —, jr inR 3 (2.4) 

Here Kq denotes the Macdonald (or modified Hankel) function. Using the boundary conditions 
in (|2.ip , we can express the kernel of interest as 

Go = G Z (x^) ~ Hx,x',^) (2-5) 



1 Here and in the following we use the convention in which the negative real axis represents the cut of the 
square root y'z in the complex plane. 
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where Q z is defined by (|2.3p - (|2.4|) and h(x,gf, yfz) solves the boundary value problem 
f (-A + z)h(x,x',^ = 

I for any x[ e Q (2.6) 

{ h(x,x',^/z)\^ m = Q z (x,x!)\^ gu 

In the next step we construct the operator which will be our main object perturbing — Aq 
by a single point interaction with the support at a point a; € f2. Such Hamiltonians are defined 
by the self-adjoint extensions of the symmetric operator 

f D(H ) = {^eH 2 n H*(n)\ Vfeo) = 0} 

Following the von Neumann theory AG63, RS75 , we observe that the restriction (|2.7[) has defi- 
ciency indices (1, 1); consequently, we arrive at a one-parameter family of self-adjoint operators 
H a . For fixed a £ 1, A £ C\R and a domain 17 C M 2 we have 



D{H a ) = |<0 e L 2 {Si) ^ = A + ^ o A (-^o) ; ^ x eH 2 nH^(fi), 

A feo) = ^ (a-lnVX- 27r% ,^, VX)) | (2.8) 

while for f2 C K 3 the operator domain is 

£>(fT a ) = |v e i 2 (r>)| : V = <f> x + qQo(,;x ), <t> x e H 2 n £#(0), 

A fe o ) = g^+^ + /i(x o ,x o ,x/A)^ | (2.9) 

In both cases the parameter A determines a representation of the operator domain, roughly 
speaking, a split between the regular and singular part; for a fixed choice of A the action of H a 
is the following 

H a iP = -Ac!) X -\qg5(;X ). (2.10) 

It is convenient to include into this description also infinite values of the parameter a, in which 
case the coefficient q of the singular part — sometimes referred to as the charge of the state — 
vanishes. It is equivalent to the absence of the point interaction: the domain reduces in this case 
to H 2 n Hq(H) and the self-adjoint extensions corresponding to a = ±oo is identified with the 
unperturbed operator — A^. 

Making use of Krein's formula [AG63J the action of the resolvent, i?" = (H a +z) _1 on L 2 (fl), 
can be expressed as a rank-one perturbation of its 'free' counterpart (—A + z) _1 , specifically 

B%<p - (-A + zT\ + Qz ((-A + (&) So(',2o) ■ (2- 11 ) 

In this formula, of course, the first term is the regular part of the function R"<f S D(H a ) 
while the value q z ((— A + z)" 1 ^) (x ) denotes the corresponding charge. Using the boundary 
conditions in (|2.8p and (|2.9p we can identify the coefficient q z with 

q z = (a-lnVz-27r/i( £o ,x , % /i)) _1 , !!cl 2 (2.12) 
+ ^+h(x ,x o ,V^)) , VcR 3 (2.13) 
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Furthermore, the Fourier expansion 

A„ + z 



(-A + z)-V= £ (2.14) 



n£N 

fc<7V„ 

yields the following explicit expression 

R * (fi = V- (^^n,fc) ^ fc | ^ (y,-0n,fc) VVt.fcfep) ^ v ^ ^ 

fc<JV„ fc<JV„ 

for Slcl 2 , and its counterpart 

2. XT^^ fc+ ^ A„ + z a+ V£ + Wx * ,/I^ o0 '- o) (2 ' 16) 

fc<JV„ fc<Af„ 

for !lcl 3 . 

3 The principal eigenvalue of H a 

As usual, the Krein formula allows us to determine the spectrum through the denominator of 
the perturbation term. In particular, it follows from the resolvent equations (|2.15p - (|2.16|) that 
the spectrum of H a is formed by the solutions of the equations 

a-lnV^;-27r/i(x , £o ,V^) = 0, SlcM 2 

(3.1) 

a+^2 + /i(x , £o ,V^) = 0, CR 3 

to which one has to add eigenvalues of — A^, the degenerate ones in any case and the non- 
degenerate ones, Xf t , provided that the corresponding eigenfunction ip n satisfies the condition 
(see e.g. in |BFM07j ) 

Vrifeo) = (3.2) 
This does not concern, however, the bottom of the spectrum which we are interested in here, 
because the ground state of — A^> is non-degenerate and can be represented by a positive function. 

The subject of this section is the principal eigenvalue of the point-interaction operator H a . 
We are going to show, in particular, that for any ael there exists a unique simple eigenvalue 
of H a below the spectral threshold of — A®. As a preliminary, we need to characterize the 
derivatives of h(x ,x , y/z) w.r.t. the variable y/z. 

Lemma 3.1 Let z be a positive real number; the function y t— > h(x,x',y), defined by \2. 6\) with 
y := y/z, satisfies the conditions 



n c 



■/■>> 



l(l + 2ird v h(x!,x!,y^ >0, 

$(& + dyh(x!,x!,v))>o, nc 

Furthermore, for any z < and h(x,x' ,iy), defined by k2. 6}) with y := yjzf, we have 

| + 2tt Imh(x!,x!,iy) = 0, SlcM 2 

jL+hnh{x',x',iy) = 0, fid 3 



(3.3) 



(3.4) 
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and 



A + 2nd y Reh(x',x',iy) < 0, SlcM 2 

d y Reh(x!,x!,iy) < 0, Slcl 3 

Proof. We start with the 3D case. Let Zj, j = 1, 2, be a pair of positive values; setting yj 
we get from equation (|2.5|) in combination with the first resolvent formula the relation 



(3.5) 



{vl-vl) (Gf,Gf 



= lim 

L 2 (Q) x-^x' 



gyi(x,x!)-g y i(x,x!) + h(x',^, yi )-h(x!,x!,y 2 ), (3.6) 



where (•, -)L 2 (m denotes the scalar product in L 2 {Q). The limit at the r.h.s. is easily seen to be 



lim 

x — >x' 



Q V2 (x,x')-gy^{x,x') = — ( yi -y 2 ) 
J 47T 



(3.7) 



Substituting this expression into (|3.6|) we get 



_ J_ h(x',x',yi) - h{x',x',y 2 ) 
L 2 (n) 4tt 



yi - 2/2 



and consequently, in the limit y\ — > y 2 we arrive at 

1 



4-7T 



+ d y h(x',x',y 2 ) = 2y 2 Q y J 



L 2 (Q.) 



> 0. 



On the other hand, for negative values of z a similar argument yields 



(3.8) 



(3.9) 



\ /L 2 (fi) x-^x' 



Q yi (x,x)-[G V2 j (x,x) -h(x,x,iyi)+h*(x,x,iy 2 ) 

(3.10) 



with the asterisk denoting complex conjugation and yj := u\zj\. The limit at the r.h.s. is 



lim 



g-yi(x,x')-(g-y*) (x,x') = -—{ yi + y2 ) 



(3.11) 



from which we get a relation replacing (|3.8p . namely 



I 2 2 

(2/1 - y 2 



L 2 (fi) 



(2/1+2/2) - h(g,tf,iyi) + h*{x',x',iy 2 ) . (3.12) 



The second one of the relations (|3.4p then follows from here by setting y\ — y 2 . As for the real 
part of (|3.10p given by 



\ {y\ vl) 



= lim 



Ga Vl ,Qo V *) +(Go y \Go Vl 

U U /L 2 (S1) V U U /L 2 (fl)_ 

Re (a, x') - Re (x, x')] - Re h(x',x', i yi) + Re x', i y 2 ) (3.13) 



we notice that the first term at the r.h.s. is in fact zero; dividing the remaining ones by y\ — y 2 
and taking the limit as y\ — ► y 2 , we arrive at 



d y Reh(x' ,x' ,iy 2 ) = 2y 2 



-y-> 



L 2 (n) 



> 0. 



(3.14) 



In the 2D case the validity of relations (|3.3|) , (|3.4p and (|3.5p can be checked following the 
same idea. Taking into account the logarithmic singularity of G z (x,x') as x — > a;' — see, e.g., 
[XS72] — we find for z > 



lim 

x — >x' 



2tt 



(In yi - In y 2 ) , 



while for z < we have 



lim 

a; — 



g- y *(x,x!) - (Q~ vl )* {x,d)\ = (-iTr + ]nyi-lny 2 ) 



(3.15) 



(3.16) 



this concludes the proof. ■ 

In the next two lemmata, we deal with the solutions of equations (|3.ip below the spectrum 



of — Aq. We are going to show that for a fixed real a there is a unique such solution. It is 
convenient to treat the 2D and 3D cases separately. 

Lemma 3.2 Let Xq denote the first eigenvalue of — Aq corresponding to the domain C R 3 . 
For any at K, the equation 



a + — h/ifeo.^o, v-£) = °! £e(-oo,A ) 

admits a unique solution, denoted £ (a), swc/i i/iai 

lim £(a) = — oo , £(— h(xn,Xn, 0)) = , 

a — > — oo 

and 

lim £(a) = Aq . 

a — >+oo 



(3.17) 



(3.18) 



(3.19) 



Proof. In order to study solutions of (|3.17p . we need to find the dependence of h(x ,x Q , \/ — £) 
on the variable \J — £. We start by considering the case £ < 0. Setting y = the equation 

(|3.17p assumes the form 

(3.20) 



— = -a-h{x ,x ti ,y) = Q, y>0 



and its solutions can be geometrically interpreted as the abscissas of the intersection points 
between the curves at the left and the right-hand side of (|3.20p . We will show that for any fixed 
choice of x S f2, h{XQ,x 10 ,y) * s a positive and strictly decreasing function of the variable y, such 
that 

lim hixnjXjQ, y) = . (3-21) 



Let us consider the boundary value problem 



(-A + y 2 ) h(x,x o ,y) = 
h{x,XQ,y)\^ dn — j-r— I 



; x E £1 



(3.22) 



xeesi 



The solution of (|3.22p is infinitely smooth in the open set Q, continuous and positive on its 
boundary. The strong maximum principle — cf. |Ev98 — in this case allows us to claim that h 
is strictly positive in f2 reaching its maximum on the boundary, 



< h(x,x ,y) < sup 



-y\x_-x \ 



for Vx e fl, y > . 



(3.23) 
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Furthermore, the derivative d y h satisfies the equation 

f (-A + y 2 ) d y h(x,x a ,y) = -2yh(x,Xo,y) 



d y h(x,x Q ,y)\ 



g -»|*-£ | 



xedQ. 



4tt 



(3.24) 



xedn 



the solution of which belongs to C 2 (ft) fl C(ft) in view of the regularity of the source term and 
the boundary value. The maximum principle — see the version given in |Br83[ Thm IX. 27] - 
in this case implies 

d y h(x,x ,y) < for Vx <= ft, y > 0. (3.25) 

In particular, the solution of (I3.24| for y = is d y h(x, x , 0) = — t=-. This characterization 
of h(x Q ,x ,y) allows us to claim that the equation (|3.20[) admits at least one solution for any 
a E (-co, — /i(xq,^o, 0)] and that the conditions (|3.18p hold. Moreover, the second one of the 
relations (|3.3[) in Lemma |3 . 1 1 implies the monotonicity of the function a <— > y(a) implicitly defined 
by (|3.20|k this grants the uniqueness of the solution. 

Next we turn to (|3.17|) for £ G (0, Aq). In this case, setting y := ^/?, the equation reads as 



a + — + ft-(x ,s ,iy) = 0. 



According to the second one of the relations (|3.4[) . this is equivalent to 

a + Re h(x , x , i y) = , 
where Re /i satisfies the boundary value problem 

(— A — y 2 ) Reh(x,x ,iy) = 



-r» 7 / \ I cos y \x — x n 

Reh{x,x ,y)\^ dn = , ' 



a; £ ft 



(3.26) 



(3.27) 



(3.28) 



xedn 



It is worthwhile to notice that h(x,x , y) is not defined for y = \f\$. In particular, one can show 
thafl 

lim eh(-,Xn,WXa - e) + ipo(x n ) ii>a(-) =0, 



(3.29) 



where ipo is the principal eigenstate of the Dirichlct Laplacian and U is any subset of ft. In view 
of the boundedness of ipq and the arbitrariness of U, this relation also implies 



lim \h{x ,x ,iy)\ 



(3.30) 



Using this result together with the conditions (|3.4p and (|3.5I) of Lemma 13.11 we conclude that 
Re h(x ,x Q , iy) is a strictly decreasing function of y S (0, \/Ao) whose behavior for y — > y/Xo is 
given by 

lim Re h(x , x , i y) — — oo (3.31) 

Summing up this discussion, the equation ()3.27p has a unique positive solution y — y(a) for any 
a S (— ^XqjXq, 0), +oo), and this solution asymptotically approaches the value Ao as a — > +oo; 
in combination with the first part this concludes the proof of the lemma. ■ 



2 Relation l|3,29| l easily follows, e.g., from Lemma 2 in [BFM07 . 
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Next we deal with the eigenvalue equation in the two-dimensional case. Recall that the free 
Green's function related to this problem is the modified Bessel function Kq, which is strictly 
positive and convex in R + and admits the following representation |AS72j 

+00 2n 

Ko(x) = - fin ~ + 7} Io(x) + J2 c « 7^77^ > ( 3 ' 32 ) 



n=l 



where 7 ~ 0.577 is the Euler-Mascheroni constant, c„ = Y^k=i an d Io( x ) is the other modified 
Bessel function given by the series 



+00 2n 

In the following proof we will make use of the asymptotic properties of Kq, 

lim K Q {x) = +00 , lim K (x) = . (3.34) 

Lemma 3.3 Let Xq denote the first eigenvalue of — Aq corresponding to the domain ft C K 2 . 
For any a£l, the equation 



a — In \J — £ — 2irh(x , x , \J — £) = , £e(— oo,A ), (3.35) 

admits an unique solution, denoted £(a), such that 

lim £(a) = -oo, £(/(xo, So , 0)) = , (3.36) 

a — >+oo 

and 

lim £(a) = A , (3.37) 

a — > — 00 

where f(x, x , \f— 5) m i/ie second one of the relations \3. 36]) is a C + °°(fl) DC '(Q) -regular func- 
tion of both the spatial variables defined by 



fix, x o ,V-0 = ^h(x, x o ,V-0+ In V -^o(V -C la- Sol), C<A . (3.38) 



Proof. The argument follows the same line as in Lemma [3721 the main difference coming from 
the specific form of the Green functions in two dimensions. For £ < and y := \/|£f, equation 
(|3.35p reads 

a = lny + 2irh(x ,x Q ,y) , (3.39) 

where h(x,x , -y/y), the solution of the boundary value problem 
(-A + y 2 ) h(x,x o ,y) = 



Ms,So>J/)Le9f2 = ^ K o(y ls-Sol)L 



; £ e n (3.40) 



is strictly positive in f2, as it follows from the maximum principle and the positivity of the 
boundary values. Consequently, the r.h.s. of (|3 . 39[) diverges as y — > +00. Moreover, from (|3.3p 
we know that lny + 2ir h(x , x , y) is strictly increasing as a function of y in the whole K + . In 
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order to study the behaviour for y — > 0, we introduce the auxiliary function f(x, x ,y), defined 
in (I3.38P which solves the boundary value problem 



(-A + y 2 ) f(x,x o ,y) = 



I ; z g ft (3.41) 



Using the definition (|3.33[) . it is easy to verify that the quantity lny + 2it h(x , x , y) coincides 
with the value of / in x , thus for y — > + we have 

lim (lay + 2ir h(x a ,x ,y)) = lim f(x ,x ,y) . 

In this limit, the problem (|3 .41 [) simplifies to 
-Af(x,x o ,0) = 



; x en (3.42) 

xedn 



where the representation (|3.32p has been considered; the regularity of the boundary condition in 
(|3.42j) implies f(x,x ,Q) G C + °°(ft) n C(f2). So far we have shown that (In?/ + 27T /i(x , z , jy)) 
is a strictly increasing function of y G K + such that 



and 



lim (In y + 2ir h(x Q , x , y)) = +00 (3.43) 

y — *+oo 

lim . (\ny + 2irh(x ,x Ji) ,y)) = /fe >^o> )- ( 3 - 44 ) 

y^0+ 

From here we can conclude that the equation ()3 .35(1 admits an unique solution, £(a) < 0, for any 
a € [f(x ,x , 0), +00), which satisfies the conditions (|3.36p . 

In the case of a positive eigenvalue £ G (0, Ao), setting y = -y/f we rewrite equation (|3.35|) as 

a = In i y + 2tt h(x Q , x , i y) (3.45) 

which due to (|3.4|) is equivalent to 



a = In j/ + 2ixYleh(xQ,Xj Q ,iy) . (3.46) 

As in the 3D case, we notice that the r.h.s of ()3.46p is a strictly decreasing function, cf. (|3.5[l . 
diverging as y — > VAo; the sought conclusion easily follows. ■ 

Since the spectrum of H a is determined by the solutions of the equations ()3.ip , the above lem- 
mata have the following implication which means that in a sense point interactions in dimension 
two and three can be always regarded as "attractive" . 

Corollary 3.4 For any real a, the operator H a defined by i2.8] )- ^.10\) has a unique spectral 
point below the spectral threshold of— A^. 



4 Dependence of the principal eigenvalue on the position 
of the interaction 

Now we pass to our main topic. We will characterize the behaviour of the principal eigenvalue 
of the point-interaction Hamiltonians H a for a fixed bounded domain £1 C M d , d — 2,3, as the 
interaction site moves towards the boundary of Q. We will restrict our attention to domains 
having an interior reflection property w.r.t. a suitable hyperplane, in the following sense: 
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Definition 4.1 Consider a hyperplane P of dimension d—1 in M. d and denote by S p the mirror 
image of a set S C M. d w.r.t. P provided S P — . The domain is said to have the interior 
reflection property w.r.t. P if P n f2 and there is an open connected component fl s C Q\P 
such that fif is a proper subset of£l\Cl s . We call fi s the smaller side of ft and P an interior 
reflection hyperplane. 

To prove our main result, we need following auxiliary statement. 

Lemma 4.2 Let Gq(x,x') be defined by l£.5|) -i TO]) and z£l. For values of z above —\q, the 
following implications hold, 

z e [0,+oc) =>■ Go(x,x!) > in Q (4.1) 

and 

ze(-\ ,0):^:Reg$(x,x!)>0 in Q . (4.2) 
Proof. For C M. 3 and z € [0, +oo), the Green function 

G z o(2L,x') = - ^~ % - h{x,x>, y/z) (4.3) 

4:TT\X — X'\ 

is certainly positive in a small enough open neighbourhood B x i of the point x' due to the 
boundedness of h(x,x', \fz). Moreover, it solves the boundary value problem 

(-A + z)g$(x,x!) = inQ\B^ 

z o\ou = 0, Go\ 9B ^,>0 ! ( 4 - 4 ) 

It follows from the maximum principle that Go is strictly positive in the whole O. 

To prove the other implication in the 3D case, notice that for z € (— Ao, 0) we have 

ReG^x') = °° S "ff ~ ~ 1 - Re h(x, x', t VR) ■ (4-5) 

47T \X — X' | 

Once more we can find a suitable open neighbourhood of the point x', which we call B x > 7 where 
this function is positive. In Q\B X > Re Gq (x, x') solves the boundary value problem 



GoU = 0, GoU>0 ^ (4.6) 



{-A + z)ReGo(x,x') = in fl\B x > 

Under the condition z > — Ao we can still apply the maximum principle obtaining in this way 
Re Gq (x, x') > in O. Finally, in the 2D case the proof follows the same line with the replacement 
(|4.3p and (|4.5p by the corresponding 2D Green's function 

G^(x,x!) =K (^\x-x!\)-h(x,x!,yfi) (4.7) 

and taking into account the asymptotic properties of Kq(x) as x — > 0. ■ 

Now we are in position to prove our main result. The next theorem shows that, under the 
interior reflection conditions imposed on the domain Q, the principal eigenvalue of H a increases 
as the interaction site moves towards the boundary of the smaller side f2 s of Q. 

Theorem 4.3 Let P be an interior reflection hyperplane for the domain and denote by n the 
normal vector to P pointing towards f2 s . Assume that x £ f2 n (dfl s D P); then the principal 
eigenvalue £ of the point-interaction H a with the perturbation placed at x satisfies the condition 

n ■ V, o £ > . (4.8) 
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Proof. Consider first the 3D case. To analyze the dependence of £ on the interaction po- 
sition x , we have to distinguish between the negative and positive spectral points. If a S 
(—00, — h(x , x , 0)], then £ = £(x ) is by Lemma f3.2l a negative solution of (|3.20p . Replacing 
with y(x ) and taking the gradient w.r.t. a; in (|3.20j) we find 



(4.9) 



V^aj) ( — + d y h ) = -V^Hx^XQ^y) 



Next we consider the term ^x ( h(x ,x ,y) at the r.h.s. of the last equation; under our interior 
reflection assumptions we will show that this vector is oriented towards the smaller side of Q. 
To this aim we notice that, in view of the relations (|2.2p - (|2.4p . h(x,x',y) can be written as 



h(x,x',y) = 



-y\x—x 



4tt\x — x 



k<N n 



- 0n,fc(£ , )'0«,fc(a) 

K + y 2 



(4.10) 



for any x ^ x' . From the symmetry of this expression and the regularity h(x,x',y) it follows 
that 

Vx>h(x,x!,y)\ x=x , =x =S/ S _h(x,x',y)\^,^ (4.11) 



and 



V £o /i(£o.2o.!/) = 2 VJi(x,x!,y)\ x=x , =x 



(4.12) 



To analyze the orientation of this vector, we introduce the function u defined on the smaller part 
of Q by 

u{x,X£,y) = h(x,X£,y)-h(x p ,x Xj ,y), xett s , (4.13) 
where x p denotes the mirror image of x S tt s through the plane P. The following equation holds 

(-A + y 2 ) u = in Q s 



u lpnn = ' u \ 



dn s ndn 



e -«|x- a3 | 
4-7T \X— Xr. 



H^. P ,x ,y) 



; SjeSlnP (4.14) 



It is worthwhile to notice that the boundary value on dQ s n dfl can be identified with the value 

2 p 
of Qq (x,Xq) on the set (dQs H dfl) , indeed we have 



e -y\x-Xo\ 



(dnsndn)' 



4tt\x~x 



h(x,x ,y) 



.0 



(dnsndn) 1 



e -y\x.-xo\ 
47T \x — x n 



h(x p ,x ,y) 



dnsndn 
(4-15) 

Then it follows from (|4.ip that m is positive on dQs n c?il and by the maximum principle, u > 
holds in f2 s . In particular, u reaches its minimum on the points of the open surface P (~1 f2; the 
Hopf boundary-point lemma in this case implies 



n ■ Vx_u > for ViePflO. 
Due to the definition (|4.13p . in combination with the relation 

n ■ Vx_h(x p , x , y) ~ ~n ■ V x ph(x F , x , y) in Q s U (P D f2) 
the last inequality also implies 

2 n ■ Vxh{x, x 0l y) > for Vx eFdSl. 



(4.16) 
(4.17) 
(4.18) 
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Substituting (|4.12[) and (14.18|) into the r.h.s. of (|4.9[) and taking the projection in the direction 
of the vector n we get 

&• V^J/Cxq) [j^ + dyh) =~Ik-Vx o h{x ,x o ,y) <0. (4.19) 

The term V^ o y(a; ) at the l.h.s. of (|4.19p is related to V £o £ by 

V, Q y = --V, o e (4.20) 

-0 y -0 



from which it follows that 



1 y{i +d « h ) W ^ > ° (4 ' 21) 



y 

The sought inequality (|4.8|) follows easily from (|4.2ip taking into account the condition 

In the opposite case, a > —h(x o ,x ,0), the first spectral point of H a is a strictly positive 
solution of the equation 

a + Reh(x ,x ,i y/l) = (4.22) 

with £ < Ao — cf. (|3.27[) in Lemma [3~2l Replacing with y(x ) and taking the gradient w.r.t. 
x Q in (j4"2"2"]) . we get 

Vx a yd y Reh(x ,x ,iy) = -V £o Re h(x a , x , i y) . (4.23) 

In order to check the orientation of the vector at the r.h.s of this expression, we notice again 
that due to the symmetry of the function 

Re x!, iy)= _ £ ^fe^-fe) (4 24) 

l_ _ 1 nGN y 

k<N n 

the gradient Re/i(^ ,s ,iy) can be expressed as 

Vx Reh(x^,^,iy) = 2VxReh(x 1 x',y)\ 2 , =sil=sLo . (4.25) 
Then we follow the line of the first part of the proof introducing the function u, 

u(x,^,y)=Reh(x,x ,iy) — Reh(x p ,x JO ,iy), x<eQ s . (4.26) 
Proceeding as before and taking into account the implication (|4.2p we find easily 

n ■ V^u > (4.27) 

and 

2n- VxReh(x,x ,y) > for ViePnO; (4.28) 
substituting (|4.25|) and (|4.28|) into the r.h.s. of lj4.23|t . we conclude that 

n-^x a yd y Reh{xj ,x JQ ,iy)<Q (4.29) 

The claim ()4.8[) is then obtained from (14.29(1 by taking into account the relation 

V Ko y = -V^e (4.30) 

-0 y -0 

and the inequality (|3.5p from Lemma 13.11 This concludes the argument in the three-dimensional 
case; the two-dimensional one can be dealt with in the same way, step by step. ■ 
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5 Optimization of £(azo) 



By Theorem 14.31 the spectral threshold of the operator H a increases as the interaction position 
x moves towards the boundary of the domain £1. This result provides us with some insights on 
how to place the point-interaction centre to minimize the principal eigenvalue of the Hamiltonian 
H a . For the sake of simplicity, we begin with the case of a convex Vt. Let II be the set of all the 
hyperplanes P of interior reflection for f2; we denote by f2 Sj p the smaller part related to P G II, 
provided it exists, and by E the union 

E = |J Vs,p (5.1) 
Pen 

The following claim is a straightforward consequence of the Theorem 14.31 

Corollary 5.1 Let be an open convex domain in K d , d = 2,3, and let H a be a point-interaction 
operator in fl with the perturbation placed at x . The principal eigenvalue of H a , considered as 
a function of the interaction centre, takes its minimum value when x belongs to the open set 
Q\E. 

Proof. Notice first that the continuity of h(x,x Q ,y) implies the continuity of solutions of the 
eigenvalue equations (|3.1[) . thus the principal eigenvalue £,(x ) has at least one minimum point 
Xq 1 £ f2. We use reduction ad absurdum: assume Xq £ fl s ,p for a suitable hyperplane P G II. 
Due to the convexity of the domain, it exists another hyperplane P' £ Id parallel to P and such 
that 

x£ g dn SiP , n p' , (5.2) 

however, under this assumption Theorem 14.31 implies the inequality 

n • V £o £ > (5.3) 

from which a contradiction follows easily. ■ 

In the case of highly symmetric domains such as the interior of a circle or an ellipse in the 
plane, and similarly a ball and an interior of an ellipsoid in three dimensions, it is easy to identify 
the set E with the center of such a domain. More generally, the convexity of fl ensures the validity 
of the interior reflection property needed in Theorem 14 . 31 with respect to some hyperplane passing 
through a point sufficiently close to the boundary. In this situation the above result can be used 
to localize the optimal position of x^ in a 'central' subset of the domain fl. 

A slight generalization of the above argument lead us to an analogous rule to localize the 
minimum points of the principal eigenvalue for a point interaction within non-convex domains. 
Let P £ LT and consider the one parameter family of hyperplanes Pt , 

P t = {x + n t : x £ P, t £ [0, T]} 

(5.4) 

T = max {t £M+ : P t nn ^ 0} 

where n denotes the unit normal to P directed towards the smaller part £l s of il. We denote as 
IT the subset formed by all those hyperplanes P G II such that 

P t £ LT for V< G [Q,T] . (5.5) 

It is important to notice that in the non-convex case, to any hyperplane of interior reflection 
there may correspond more than one smaller part. Next we denote by 6 s ,p the union of all the 
smaller parts related to P, and by E' the set 

E' = (J 9 s ,p . (5.6) 
Pen' 
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Corollary 5.2 Assume that is an open domain in R d , d — 2,3, and H a is a point-interaction 
operator in Q with the perturbation placed at x Q . The principal eigenvalue of H a , regarded as a 
function of x , takes its minimum value when Xq belongs to the open set f2\E'. 

Proof. The argument is an easy modification of the proof of Corollary 1 5. II ■ 

Non-convex domains with a reasonably simple boundary such as, for instance, the union of 
two intersecting disks or a dog-bone profile in two dimensions, can be easily analyzed using 
Corollary 15. 21 It is also worthwhile to stress that the results of this section do not depend on the 
parameter a, hence the optimal placement of the point interaction with respect to the minimum 
of the principal eigenvalues can be the same irrespective of the interaction "strength" . 

Let us finally comment on he relation to the work HKK01 mentioned in the introduction. 
We have said that for a hard-wall obstacle the principal eigenvalue decreases as it moves towards 
the boundary. The difference of the two effects can be traced back to the different boundary 
conditions which characterize the operator domains in the two cases. While the hard obstacle 
is characterized by Dirichlet boundary condition, the point-interaction operator H a considered 
here can be obtained as the norm-resolvent limit of a family of sphere interactions Hamiltonians 
H a {r) with the boundary condition of a mixed type as the radius r — ► 0. In the three-dimensional 
case, for instance, the operator H a (r) is explicitly given by 

j H a (r) = -A on Q\S r 

[ D(H a (r)) = {^e H 1 (M. 3 )\ (S„V0+ " M)- = i^kp^} 

where S r denotes the sphere of radius r centered at x and {d n ip), — (d n tp)_ is the jump of 
the normal derivative of ip on the interaction surfac e |FT93[ ISh03j . Another insight into the 



difference of the two situations can be obtained from ES96 
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